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The most efficient energy sources known in the Universe are accretion disks. Those around 
black holes convert 5-40 per cent of rest-mass energy to radiation. Like water circling a 
drain, inflowing mass must lose angular momentum, presumably by vigorous turbulence in 
disks, which are essentially inviscid The origin of the turbulence is unclear. Hot disks of 
electrically conducting plasma can become turbulent by way of the linear magnetorotational 
instability Q]. Cool disks, such as the planet-forming disks of protostars, may be too poorly 
ionized for the magnetorotational instability to occur, hence essentially unmagnetized and 
linearly stable. Nonlinear hydrodynamic instability often occurs in linearly stable flows (for 
example, pipe flows) at sufficiently large Reynolds numbers. Although planet-forming disks 
have extreme Reynolds numbers, Keplerian rotation enhances their linear hydrodynamic 
stability, so the question of whether they^can be turbulent and thereby transport angular 
momentum effectively is controversial 3|, 
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151 ] . Here we 



report a laboratory experiment, demonstrating that non- magnetic quasi-Keplerian flows at 
Reynolds numbers up to millions are essentially steady. Scaled to accretion disks, rates of 
angular momentum transport lie far below astrophysical requirements. By ruling out purely 
hydrodynamic turbulence, our results indirectly support the magnetorotational instability 
as the likely cause of turbulence, even in cool disks. 

Our experiments involved a novel Taylor- Couette apparatus [if]]. The rotating liquid 
(water or a water /glycerol mixture) is confined between two concentric cylinders of radii 
Ti,T2 (ti > fi) and height h. The angular velocity of the fluid is controlled by that of the 
cylinders, Q\ and ^2- An infinitely long, steady, Taylor- Couette flow rotates as 

0(r) = a+^, (1) 

where a = (f^r 2 — ^irf)/(r| — r i) an d b = rfr^^i — (^/(r 2 . — rf). Reynolds number here 
can be defined as fir^ — ri)(Oi — Qq)/v, where v is viscosity and f = {t\ + r?)/2. The 
rotation profile (equation (TO) ensures a radially constant viscous torque —2ir puhr 3 dQ / 'dr 
for constant mass density p and constant v. Astrophysical disks are mostly Keplerian, 
meaning oc r~ 3//2 , so that decreases radially outward (d\Q\/dr < 0) while the specific 
angular momentum, |r 2 f2|, increases radially (d\r 2 Q\/dr > 0). We apply the term "quasi- 
Keplerian" to any flow satisfying these conditions, which are crucial for both hydrodynamic 
and magnetohydrodynamic linear stability j^. 

Real flows have finite length. Disks have nearly stress-free vertical boundaries, but viscous 
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stress at the vertical endcaps of laboratory flows drives secondary circulation. This may 
cause the rotation profile to deviate significantly from equation (pQ), and may even provoke 
turbulence Our apparatus incorporates two rings at each end that 

are driven independently of the cylinders (Fig{TJ). Thus we have four controllable angular 
velocities. When we choose these appropriately, secondary circulation is minimized, and 
ideal Couette profiles are closely approximated throughout the flow except within ~ 1 cm 



of the endcaps [21]. 

Most past work has taken the outer cylinder at rest (f^ = 0), so that both |^| and |r 2 f2| 



decrease radially outward. Such flows are axisymmetrically linearly unstable [la] at modest 
Reynolds number, Re. Very few experiments have studied the linearly stable regime where 
d\r 2 Q\/dr > 0, as occurs in disks. Among these few are two classic experiments of the 
1930's: in one of these, the inner cylinder was at rest (fii = 0)[l8], while in the other 
< Q1/Q2 < 1- Enhanced torques between the cylinders and other evidence of turbulence 
were reported at sufficiently large Re. These results have been cited in support of the hy- 

nn 

pothesis of nonlinear hydrodynamic turbulent transport in disks ]3, 6], notwithstanding that 
the direction of turbulent angular momentum transport, which always follows —dfl/dr on 
energetic grounds, differs in sign between these experiments and astrophysical disks. The 
so-called f3 prescription |(| derived from the above experiments has been used to model or 



interpret astronomical observations[e.g. 



221 ] . To our knowledge, the only published experi- 



ments with quasi- Kepler ian flow at relevant Reynolds numbers are those of Richard] 
and Beckley 23]. In Richard's work, transition to turbulence via nonlinear instabilities was 
studied qualitatively by a flow- visualization method. No direct measurements of angular 
momentum transport were performed. Beckley did measure torques roughly consistent with 
the (3 prescription but attributed them to secondary circulation, which was strong because 
his endcaps co-rotated with the outer cylinder and h/(r 2 — ri) was only ~ 2. 

Experimental flows studied by ourselves and others are summarized in Fig. [2j Our 
Reynolds numbers are up to 20 times larger than those previously achieved by Richard^]. 
A Laser Doppler Velocimetry (LDV) system was employed to sample the azimuthal velocity 
vg at various radial and axial locations. Both mean values, Ug = rQ, and fluctuations, 
v' e = vq — vg, were obtained. At our Reynolds numbers, linearly-unstable flows are always 
turbulent, with fluctuation levels (v' e 2 ) l l 2 jvg = 5 — 10%. They are largely insensitive to the 
endcap speeds. In contrast, quasi-Keplerian and other linearly-stable flows are sensitive to 
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the endcap boundary conditions. When the endcap speeds are adjusted to best approximate 
ideal Couette flow, the fluctuations are < 1 — 2% and indistinguishable from those of our 
solid-body flows, which are expected to be steady or laminar due to the lack of shear to 
drive turbulence. 

One hallmark of nonlinear instability is hysteresis: the transition from laminar flow to 
turbulence occurs at higher Reynolds numbers than the reverse|7|. Quiescent flows were 
gradually brought into linearly- unstable regimes by raising Qi or lowering Q 2 , then returned 
to linearly stable regimes. Significant fluctuations were found only in the linearly unstable 
regime; no hysteresis was detected. 

In the absence of magnetic fields, turbulent angular momentum transport requires cor- 



related velocity fluctuations. The radial angular momentum flux is prv' g v' r where v' r is the 
fluctuation in radial velocity, and the turbulent viscosity f tU rb is defined by equating this to 



—pu turh r 2 dQ/dr where z/ turb = (3 \r 3 dQ/dr |, so that (3 = v' g v' r / (r 2 dQ/dr) 2 is dimensionless. 
Conveniently, in a turbulent but statistically steady state with profile (P), both (3 and z/ tur b 
are radially constant. A value of (3 — (1 — 2) x 10 -5 has been deduced 6] from experiments 
lillS! with < fij/fta < I- 

We have measured the Reynolds stress directly using a synchronized dual-LDV system. 
Both v'q and v' r appear to follow Gaussian statistics (E.S. et al., manuscript in preparation), 
and random errors were reduced by averaging 10 3 to 10 4 samples. Systematic errors were 
gauged by comparison with solid body flows (Oi = Q 2 = ^3 = ^4)- Figure [3] shows 
that (3 differs indistinguishably between quasi-Keplerian and solid-body rotation, falling 
far below the proposed range [6(. A large outward Reynolds stress is detected in linearly- 
unstable flows, (3 > 10~ 3 , as expected (M.B. et al., manuscript in preparation). Even for 
linearly-stable flows, when the speeds (f^f^) of the endcaps were not adjusted properly 
to produce the ideal-Couette profile ([1]), (3 was found to be almost 10~ 4 (Fig. H]). This 
again indicates that the axial boundaries can profoundly influence linearly-stable flows, as 
previously suggested [15|, 1 19f] - 

A final point should be made about the critical Reynolds number for transition, -Re cr it, 



13, 



181 ]. the scaling 



versus gap width. Based on the experiments of Wendt and Taylor 
i?e crit — 6 x 10 5 • [(r 2 — ri)/r] 2 has been proposed^, 0, [f3, Q]. It is unclear whether this 
scaling applies to quasi-Keplerian flow since it was derived from data for f^/^i > 1. In any 
case, up to Re = 2 x 10 6 , which is about three times the proposed Re crit since (r 2 — r 1 )/r ~ 1 
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in our device, we have seen no signs of rising fluctuation levels or Reynolds stress (Fig. H]). 

Therefore, we have shown that purely hydrodynamic quasi-Keplerian flows, under proper 
boundary conditions and at large enough Reynolds numbers, cannot transport angular mo- 
mentum at astrophysically relevant rates. 

Of course, it could be argued that our maximum Re, which only barely exceeds some 
theoretical estimates of -Re cr i t [3] , is still not large enough for transition. Or it could be that 
transition has occurred, but that the transport is too small for us to detect. To extrapolate 
from Re < 2xl0 6 to a typical astrophysical value > 10 12 , we rely on the empirical observation 
that for Re > Re cr it, the Reynolds number based on the turbulent rather than the molecular 
viscosity is approximately independent of Re itself [24J ]. It follows that u turh « LU/Re CT i t 
for Re > Re CT it, where L and U are the characteristic size and velocity of the flow. It is 
common knowledge among civil engineers that this is true of flow in pipes, for example. If 
it is true of rotating shear flow, as theoretical arguments suggest it should be , then (3 at 
Re > 10 12 should be comparable to what we find at Re < 2 x 10 6 , namely (3 < 6.2 x 1CT 6 (at 
2 s.d., or 98% confidence; see FigfJ] caption) — whether or not we have crossed the threshold 
of transition. 

Lastly, it is useful to relate the above upper bound for /3 to the more commonly used 
Shakura-Sunyaev a parameter [lj]: z/ tur b = aVth? . We replace this by z/ tur b = c^( r 2 — r i) 2 
since r 2 — r\ is smaller than h in our experiment, and we presume that the dominant 
turbulent eddies scale with the smallest dimension of the flow (h <C r in most disks). Then 
a = j3qr 2 / (r 2 — ri) 2 ~ (3q (see FigfS] caption for a definition of q), and thus our results imply 
a similar upper bound for a as for (3 in purely hydrodynamic disks, whereas protostellar- 
disk lifetimes and accretion rates indicate[25| a > 1CT 3 . Although some have suggested 
that complications such as vertical or radial stratification may yet lead to essentially linear 
nonaxisymmetric hydrodynamic instabilities 26j, |27|, our belief is that such nonaxisymmetric 
linear instabilities depend upon radial boundaries and hence are not generally important in 



thin disks 



28 



29(. If this is correct, then by default, the magnetorotational instability 



appears to be the only plausible source of accretion disk turbulence. 
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FIG. 1: Experimental setup. A rotating fluid (water or a water /glycerol mixture) of height 
h = 27.86cm is confined between two concentric cylinders of radii n = 7.06cm and T2 = 20.30cm, 
which rotate at rates of f2i and O2, respectively. Two novel features distinguish this apparatus from 
conventional Taylor-Couette experiments. First, secondary circulation is controlled by dividing 
each endcap into two independently driven rings. Opposing rings at top and bottom are driven at 
the same selectable angular velocity Q3 (inner rings) and (outer rings). Traditionally, a large 
aspect ratio T = /i/( r 2 — r i) is used to reduce the secondary circulation. However, at T ~ 25 by 



Richard [7] and even at T ~ 100 by Taylor 



la ], end effects were reported to be significant when the 



endcaps co-rotated with one of the cylinders. Even when the endcaps were divided into two rings, 
but with each affixed to one cylinder , residual secondary circulation may have facilitated 

observed turbulent transitions [I5I When and ^4 are appropriately chosen, secondary 

circulation is minimized and ideal Couette profiles are well approximated. 2l|]. A second novel 



feature is access to rotation profiles on both sides of marginal linear stability at Reynolds numbers 

as large as 10 6 [also see Fig.(j2])]. When the specific angular momentum, r 2 ft, decreases with 

increasing r, the Rayleigh stability criterion 3(| is violated, and thus the flow is linearly unstable 

when Reynolds number exceeds a critical value When d\r 2 Q\/dr > but d\Q\/dr < (as in 

disks, where oc r" ~ 3 / 2 ), the flow is quasi-Keplerian and known to be linearly axisymmetrically 

stable. All major components of the apparatus were precisely machined and balanced, and except 

for the inner cylinder and rotating shafts, are made of clear acrylic to facilitate visual and laser 
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diagnostics. 
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FIG. 2: Experimentally studied Taylor-Couette flows. Axes are Reynolds numbers based 
on the inner and outer cylinders: Re\^ = ^i,2 r i,2( r 2 — r i)/v. Asterisks mark Rayleigh-unstable 
flows; squares, quasi-Keplerian flows, i.e. d\£l\/dr < but d\r 2 Q,\/dr > 0; diamonds, solid-body 
flows (dil/dr = 0]; crosses for the inner cylinder at rest; triangles for flows explored in previous 



experiments 
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19(; and the rectangular box for the parameter regime explored by Richard 7j]. 



Dashed lines denote constant values of q = —dlnQ/dlnr at r = 17 cm, where most of our 
measurements of Reynolds stress were performed (Figj3]). Rayleigh-unstable flows exhibit 5 — 10% 
fluctuations that are insensitive to the end-ring speeds. Quasi-keplerian and other Rayleigh-stable 
flows are more sensitive. For example, when the end rings are fixed to the cylinders, fluctuations 
up to 4 — 8% occur. When the end-ring speeds are adjusted so that the ideal-Couette profile is 
restored, fluctuations are < 1 — 2% and indistinguishable from those of our solid-body flows, except 
within a few cm of the boundary. Reducing Re by a factor ~ 18, using an admixture of glycerol, 
actually increases the fluctuation level for the same (nearly ideal-Couette) profile. We interpret 
this to mean that the residual unsteady secondary circulation penetrates deeper into the bulk flow 
at lower Re. We infer from this that the experiments by Richard[?J may have been affected by the 
endcaps, although his ratio of height to gap width exceeded ours. His endcaps were split but fixed 
to the cylinders. 
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□ [400,146,40,53]rpm 
* [53,53,53,53]rpm 




FIG. 3: Experimentally measured Reynolds stress versus height in a quasi-Keplerian 
profile. Here f3 = v d v' r /(vg 2 q' 2 )sign.(qQ), where q = — dlnQ/dlnr: q = +3/2 in Keplerian disks, 
1-2 < q < 1.9 in our quasi-Keplerian flows. Square symbols connected by solid line were taken 
at q = 1.86 (see FigJ2J. Starred points connected by dotted line are a solid-body case (q = 0), 

which should be nonturbulent and therefore serves as a control for systematic errors. The measured 

J 

values fall far below the range of [3 proposed by Richard & Zahn 6], shown with horizontal dotted 
lines. The measurements were performed using a synchronized, dual laser-Doppler-velocimetry 
(LDV) system, which allows simultaneous detection of both components of vg and v r . The laser 
beams enter the fluid vertically through the acrylic endcaps from below (see Fig{TJ), and Z is the 
height above lower endcap. Rotation speeds [fii, Q3, fi4, f^] are shown, where O3 and are the 
angular velocities of the inner and outer end rings. Apparently gaussian deviations amounting 
to ~ (1 — 2%) of v$ were observed from sample to sample, representing measurement errors and 
perhaps true fluctuations. Each data point is the mean of f3 computed from 3, 000 — 10, 000 samples. 
Error bars represent the standard deviations. 
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FIG. 4: Dimensionless Reynolds stress at Reynolds numbers up to 2 x 10 . At the 

largest Re, the data show no sign that f3[= v'gv' r /(ljg 2 q 2 )] or the relative fluctuations v' 2 /(TJg) 2 
themselves increase with Re in quasi-Keplerian flows. Here q = — d In Vt/d In r is calculated from 
the mean vg(r) profiles. Systematic errors in (3 have been accounted for by reference to identical 
measurements in solid-body flows. Different sizes in error bars (showing s.d.) for the 6 quasi- 
Keplerian flows with optimum boundary conditions (squares at Re > 10 5 ) are largely due to 
different numbers of measurement samples, N. Averaging over these 6 points, weighted by vN, 
results in (3 = 0.72 x 10~ 6 with s.d. of 2.7 x 10~ 6 or (3 < 6.2 x 10~ 6 at 2 s.d. (98% confidence). 
When the end-ring speeds are not optimized to produce Couette profiles ([T]), however, j3 exhibits 
large values (asterisks). When optimal speeds are used but glycerol is added to the water to 
reduce the Reynolds number, larger (3 values are also seen (squares at Re < 10 5 ), possibly due to 
stronger residual secondary circulation. For reference, the solid line is the transport expected due 
to molecular viscosity: (3 V - 1SC = v/[f 3 (Q,2 — ^i)/(^2 — n)] = Re~ l {r2 — r\) 2 /f 2 . The proposed f3 lies 
between the horizontal dotted lines [6]. 
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